Metal enriched gaseous halos around distant radio galaxies: 
Clues to feedback in galaxy formation 1 
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ABSTRACT 

We present the results of an optical and near-IR spectroscopic study of giant nebular emission 
line halos associated with three z > 3 radio galaxies, 4C 41.17, 4C 60.07 and B2 0902+34. 
Previous deep narrow band Lya imaging had revealed complex morphologies with sizes up to 
100 kpc), possibly connected to outflows and AGN feedback from the central regions. The outer 
regions of these halos show quiet kinematics with typical velocity dispersions of a few hundred 
kms -1 , and velocity shears that can mostly be interpreted as being due to rotation. The inner 
regions show shocked cocoons of gas closely associated with the radio lobes. These display 
disturbed kinematics and have expansion velocities and/or velocity dispersions >1000kms -1 . 
The core region is chemically evolved, and we also find spectroscopic evidence for the ejection of 
enriched material in 4C 41.17 up to a distance of «60kpc along the radio-axis. The dynamical 
structures traced in the Lya line are, in most cases, closely echoed in the Carbon and Oxygen 
lines. This shows that the Lya line is produced in a highly clumped medium of small filling 
factor, and can therefore be used as a tracer of the dynamics of HzRGs. We conclude that these 
HzRGs are undergoing a final jet-induced phase of star formation with ejection of most of their 
interstellar medium before becoming "red and dead" Elliptical galaxies. 



Subject headings: galaxies: active — 
(B2 0902+34, 4C 60.07, 4C 41.17) 



;alaxies: formation — galaxies: high-redshift 
quasars: emission lines 



galaxies: individual 



1. Introduction 
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There is compelling evidence that in galax- 
ies the formation and evolution of the central 
stellar bulge and the ma ssive nuclear black hole 

19981 : 



are intimately related 



Magorrian ct al 
GebhardtJ 12000 ) 



Ferrarese fc Merrittl l2000t iGebhardtJ 12000ft . Un- 
derstanding this coevolution of galaxy spheroids 
and their central massive black holes is one of the 
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major outstanding issues in modern cosmology. 
Because they are readily located by their ultra- 
steep spectrum radio properties, the high redshift 
radio galaxies (HzRGs; z > 3) provide an efficient 
means to locate and study the environment and 
physics of newly-forming galaxies. There are at 
least two reasons why HzRGs are key in attempts 
to understand the physical processes involved. 
First, they rank among the most luminous, largest, 
and most massive galaxies known i n the early Uni- 
verse (e.g., iDe Breuck et al. 2002 ). Secondly, we 
observe them early in the epoch of galaxy for- 
mation at a time when their super-massive black 
holes (SMBHs) are highly active, and while their 
relativistic jets are interacting most strongly on 
their host galaxies. 

Within the framework of the standard A Cold 
Dark Matter (A-CDM) scenario, it is believed that 
massive galaxies grow in a hierarchical fashion 
through the merging of smaller stellar and dark 
matter halo objects. Whether their central black 
holes grow in simil ar fashion or whether they are 
primordial objects dLoeblll993t ISiik fc Reedll99i 
Kauffmann k. HaehneltJ 12000 ). the energetic out- 
flows and the ionizing radiation from these cen- 
tral active SMBHs are expected to profoundly in- 
fluence the evolution both of their parent galax- 
ies and of the surrounding environment. Recent 
models describing the formation of massive galax- 
ies and clusters provide furt her evidence for the 
importance of this feedback (Benso n et al. 2003; 



voirs enable us to study in detail the feedback 
processes occurring at high redshift. 

2003h 



In an earlier paper (jReuland et al 



ISpringel. Pi Matteo. fc Hernquistl l2005h . Addi 



tionally, outflows could provide a source for the 
chemical enrichment of the inter galactic medium 
seen at high redshifts (e.g., Rauch, Sargent, & 
Barlow 2001; Aguirre et al. 2001). While there 
is consensus that feedback in some form must be 
important, there is an ongoing debate about which 
manifestations dominate and the scales on which 
they operate. 

Because they are themselves massive galax- 
ies in the early universe, HzRGs are located 
in the regions of large matter overdensities and 
within regions rich in forming galaxies. Indeed, 
many are known to be embedded in regions of 
dense interstellar medium (ISM) and in environ- 



ters (e.g., McCarthy 


19931: Ivan Oiik et al. 
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1997: Athreva et all 


1998: Pentericci et al. 


2000 


Papadopoulos et alJ 20001). These easeous 
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described the results from narrow-band imaging 
observations of unprecedented sensitivity of three 
HzRGs (B2 0902+34, 4C 60.07, and 4C 41.17). 
These images were obtained at the Keck II 10 m 
telescope using custom-made, high-throughput in- 
terference filters with bandpasses centered at the 
redshifted Lya line. The observations revealed 
very luminous (-^Lya 2 10 45 ergs -1 ) and extended 
(«200 kpc) emission line nebulae with spectacu- 
lar features not previously seen, such as long fila- 
mentary structures, ionization cones, and multiple 
sharply bounded regions of enhanced emission, all 
indicative of strong interactions as expected from 
the scenario painted above. We argued that these 
extended Lya nebulae might represent gas cooling 
in massive CDM halos, supplying new material for 
the continued growth of the galaxies at their cen- 
ter. This feeding could then be responsible for 
feedback mechanisms through radio jets, super- 
nova explosions, and radiation pressure from the 
AGN, resulting in large scale outflows. The ex- 
tended X-ray halo around 4C 41.17 provides fur- 
ther evidence for a highly interactiv e environment 
in such systems (jScharf et al.ll2003h . 

However, with only morphological data, impor- 
tant issues concerning possible origins for the fil- 
aments and large-scale structures well beyond the 
radio sources could not be tackled. Because Lya 
is a resonance line, it is also important to re- 
solve the question of whether the nebulae repre- 
sent truly extended ionized halos, or rather, are 
due to Lya photons which are produced near the 
central AGN and are scattered off neutral hydro- 
gen halos. Many other questions concerning the 
structure, origin, and fate of the emission-line gas 
remain. For example: is the gas organized in 
shells, filaments, or cloudlets? What is the source 
of ionization and what is the chemical composition 
of the gas? Can the outflows regulate the growth 
of galactic bulge and black-hole? Can they expel 
metals from the deep potential wells to enrich the 
intergalactic and inter-cluster medium? 

In order to better understand the kinematics, 
abundances, and ionization mechanisms of these 
halos we obtained optical and near-infrared spec- 
tra with the Keck telescopes to measure the ex- 
tent, intensity, and kinematics of Lya and other 
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non-resonance lines, [He II], [O II] and [O III], at 
various position angles across the nebulae. This 
paper discusses the results of these spectroscopic 
observations. 

The structure of the paper is as follows: The 
sample selection, observations and data analysis 
are described in §2. In §3 we present results for 
individual objects. §4 is a discussion of the ob- 
served kinematics and line ratios of the halos. The 
implications for the origin and fate of the emission 
line halos are discussed in §5 and our conclusions 
are given in §6. 

Here we adopt the concordance cosmological 
parameters = 0.27, £1a = 0.73 and H = 
71 km s _1 Mpc -1 . The age of the Universe is 
1.7 - 1.9/ifj 1 Gyr at the redshifts (z = 3.4 - 3.8) 
of our galaxies, and the angular-to-linear transfor- 
mation is «7.4ft,f 1 1 kpc arcsec -1 . 

2. Observations and data analysis 
2.1. Sample Selection 

The three objects selected for the spectro- 
scopic follow-up observations are shown in Ta- 
ble [1] with their positions and adopted redshifts. 
They were selected from amongst the galaxies ob- 
served in the course of our Keck imaging program 
(jReuland et al.ll2003t ). The reasons for their inclu- 
sion in the spectroscopic program are summarized 
as follows. 

4C 41.17 at z = 3.8, was one of the first HzRGs 
to be discovered (jChambers. Milev fe van Breugel 
1990) and for ma ny purposes serv es as an archetype 
HzRG. Optical (|Dev et all Il997h as well as sub- 
mm wavelength o bservations (jDunlop et al.lll994 ; 



Ivison et al.ll2000f) have shown that it is a mas 



sive forming galaxy with a star formation rate 
of up to several thousand MQyr -1 . Recently, 
very extended X-ray emission was found around 
4C 41.17, which follows the Lya morphology 



closely (jScharf et al.ll2003l ) 

We selected 4C 60.07 (z 



3.8: Chambers et al 



19961 : iRottgering et al.l 119971 ) because it shows 



both spa tially and kinematically r esolved CO 



emiss ion (|Papadopoulos et al.l 120001: iGreve et al.. 
20051 ). Interestingly, in |ReulandeLalJ (j2003j) it 



was found that the Lya halo has a very extended 
(76/iy/ kpc) filament which appears orthogonal to 
the major axis of the CO and dust emission. 



B2 0902+34 (z = 3.4: iLinvl [19881 ) was selected 
because it is thought to be a protogalaxy, dom- 
inate d by young stars (jEisenhardt. fc Dickinson! 
1992p. So far, it is one of only a handful of 
HzRGs for which neutral hydrogen has been de- 
tected in absorption agai n st the r adio continuum 
( Uson. Baeri fc CornwelJ Il99lt ICodv fc Braunl 
2003f ). 



2.2. Optical and Near-Infrared Spectroscopy 

Figure [T] shows the different slit positions used 
in the program overlaid on contour representations 
of the narrow-band Lya images of the galaxies. In 
contrast to most previous spectroscopic studies of 
HzRGs, most of the slits were not placed directly 
along the radio axis. The specific instrumental 
setups are given in Table El The data reduction 
techniques are described below. 

2.2.1. Optical Observations: LRIS 

Most of the optical observations were carried 
out using the Low-Resolut ion Imaging Spectrome- 
ter (LRIS; |Oke et al ]|l995l) at the Cassegrain focus 
on the 10 m Keck I Telescope. The data were col- 
lected with various instrumental setups using both 
the long-slit mode and multi-slit masks designed to 
obtain spectra for sa 15 targets in the field simul- 
taneously (as part of a survey looking for associ- 
ated galaxies in the proto-cluser; Croft et al. 2006 
in preparation) . The red-sensitive LRIS-R camera 
was employed. This uses a Tektronix 2048x2048 
CCD detector with a pixel scale of 0'.'215 pixel -1 . 

All of the spectroscopic reductions were per- 
formed using st andard met hods and the NOAO 
IRAfi package jTodJ ll223) ■ Skylines were used 
to improve the first order wavelength calibration 
based on arc spectra to better than 0.3 A rms. The 
instrumental resolution was measured from the 
unblended skylines. Flux calibrations were per- 
formed using observations o f standard stars suc h 
as Feige 110 and Feige 34 (|Massev et al.lll988h . 



The extended emission of 4C 41.17 filled the nar- 
row slits of the multi-object spectroscopic pro- 
gram, rendering accurate sky subtraction difficult. 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Fig. 1. — Contour representations of the Lya emission line halos around 4C 41.17, 4C 60.07, and B2 0902+34 
(left, middle, right) with the various PAs at which they were studied overlaid. In this figure, north is at 
the top and east at the left. The solid and dashed lines correspond with slit positions of the optical and 
near-IR spectroscopy respectively and the dotted lines indicate the positions of the radio cores. The bold 
lines represent the direction of the radio axes. In the case of 4C 41.17 two axes are shown, corresponding 
to the inner and outer radio lobes. For each object the contours indica te observed surface b rightnesses of 



6.7 x 10~ 19 x (6,12,25,50,100,200,400,800) erg s" 1 cm -2 arcsec -2 ; see lReuland et all (|2003h 



This does not seriously affect the kinematics and 
relative fluxes of interest for this paper. 

2.2.2. Optical Observations: ESI 

One set of observations along the filament 
of 4C 60.07 was made during the night of UT 
2001 February 25, usin g the Echelle Spectro- 
graph and Imager (ESI; ISheinis et aL 2000t ) at 
the Cassegrain focus of the Keck II 10 m telescope 
in low-dispersion mode. The detector used is a 
high-resistivity MIT-Lincoln Labs 2048 x 4096 
CCD with a plate scale of 0'.'154 pixel -1 . Expo- 
sures were broken into integrations of 1800 seconds 
each, one of which had to be halved because of 
time constraints. We performed 6" offsets between 
each integration. The data were rotated over an 
angle depending on the position of the object on 
the slit, to align the dispersion axes. The wave- 
length calibration varies slightly with slit position, 
hence we shifted the spectra along the dispersion 
axes. This first order approximation is sufficiently 
accurate for the region of interest 4300A-9200A 
(corresponding to approximately 900A— 1920A in 
the rest frame). Subsequent data reduction was 
done using standard methods in IRAF. 



2.2.3. Near-IR Spectroscopy: NIRSPEC 

The near-infrared spectra were obtained us- 
ing the 10 m Keck II Telescope with its Near- 
Infrared Spectrograph (NIRSPEC; iMcLean et"al 



19981 ). The slit dimensions were 0.76" x 42" slits 
giving low-resolution (R w 1400—1900) spectra in 
wavelength ranges chosen to include the H/3, [O II] 
and [O III] lines of the target galaxies (see Table 
[2] for details). In this low-resolution mode, the 
1024x1024 ALADDIN InSb detector has a plate 
scale of 0'.'143 pixel -1 . We obtained sets of 900s 
integrations each with w 5 — 10" spatial offsets 
between exposures. 

The NIRSPEC spectra need to be corrected 
for the spectral curvature and spatial distortions 
caused by the high-throughput optics. A general 
correction would require rectification onto a slit 
position-wavelength grid based on a wavelength 
solution from skylines and coadded exposures of 
a standard star. However, since no continuum is 
apparent in our spectra, we have only extracted 
small regions lying close to the emission lines. This 
approach requires only a simple rotation over an 
angle which depends on the wavelength of inter- 
est to provide a local calibration of the wavelength 
along the slit. 

The data were flat-fielded and corrected for cos- 
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mic rays and bad pixels in the standard fashion. In 
order to remove the strong near-IR skylines, and a 
sky frame scaled to the brightness of unsaturated 
sky lines near the emission line of interest was sub- 
tracted. Subsequently, the frames were cropped, 
rotated, and coadded. Flux calibration was done 
with standard stars of spectral type AOV, B3, and 
G4, and was consistent to within 10%. 

2.3. Data analysis 

For the data analysis the spectra were regis- 
tered i n position with radio maps f r om the liter- 
ature dCarilli. Owen fc Harris! Il994 ICarillilll995l: 
Carilli et al.l Il997h . The zero po nts of the spa- 
tial scales correspond with the radio core. This 
was achieved by identifying the core with the cen- 
troid of the continuum emission. If no continuum 
emission was visible, we identify the core with ei- 
ther the spatial region that shows the broadest line 
emission, or by bootstrapping to spectra for which 
the core could be reliably identified. These results 
were then checked with the narrow-band images, 
correlating with the peaks and dips in the observed 
intensity. This resulted in a total uncertainty less 
than 0.5" in the relative spatial offsets. 

The velocity scale used in the analysis is rel- 
ative to the systemic velocity derived from the 
He II A1640 line (the adopted systemic redshifts 
are given in Table [l}. The kinematic informa- 
tion was obtained from the 2-D spectra using a 
program written in IDL making use of the pub- 
licly available fitting routine MPFIlll. The data 
were coadded within apertures matched to the 
seeing in order to increase the signal-to-noise ra- 
tio and ensure that the extracted spectra are not 
correlated. We then determined the velocity cen- 
troid, FWHM (corrected for instrumental broad- 
ening measured from unsaturated skylines, and as- 
suming that the widths add in quadrature). The 
peak and baseline fluxes were determined by fit- 
ting a single Gaussian and baseline to each trace. 
Single Gaussians provide a good fit to the outer 
regions and in order to allow a direct comparison 
with the central regions similarly, we have cho- 
sen to treat all regions consistently using the sim- 
pler approach. The single Gaussian decomposition 
for high surface brightness emission may break 
down when these regions are embedded in more 
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Fig. 2. — Grayscale representations of the 2-D 
spectra of 4C 41.17 centered at the Lya, [O II] 
A3727, H/3, and [O III] A5007 lines. The veloci- 
ties are indicated relative to the systemic redshift 
determined from the He II line. The zero points 
of the spatial scales correspond with the position 
of the radio core as identified with the broad line 
and continuum emission. Note that the peaks of 
the line emission c orresponding with the position 
of radio knot B2 in lCarilli. Owen fc: Harris] (|l994f) 
are 1" E of the nucleus. It can be seen (most 
easily in the [O III] A5007 spectrum) that there 
are 0.5" offsets between the broad lines and the 
central depression. The [O II] spectrum is rich 
in very strong skylines. These have been masked 
off to better show the structure of the [O II] line 
emission. 

smoothly varying large scale low surface bright- 
ness envelopes. 

3. Dynamical Results 

Figures 03 and [5] show 2-D optical and near- 
IR spectra of 4C 41.17, 4C 60.07, and B2 0902+34 
centered at the emission lines most relevant to our 
discussion. It is immediately obvious, that many 
of the lines are very broad, and show strong spatial 
variation both in their velocity centroids and their 
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Fig. 3. — Relative velocities (solid lines), velocity 
dispersion (bars), and normalized surface bright- 
ness profiles (dashed lines) as determined from 
spectra with slit positions parallel to the radio axis 
for the Lya, [O II], and [O III] emission lines of 
4C 41.17. Top to bottom: Lya along PA=76°, 
Lya for PA=81°, [O II] for PA=70°, and [O III] 
for PA=68°. The spatial zeropoints correspond 
with the position of the radio core and the dot- 
ted lines represent the position of the radio lobes 
projected on the slit. East is at the right in all 
of these panels. The bar and symbol size indicate 
the la uncertainties on the measurements. Note 
that the fits to the [O II] emission are affected by 
the strong skylines and that the slit along PA=81 ° 
did not go through the radio core. 



FWHMs. In some cases continuum emission was 
also detected, but only from the central regions, 
within 3 "of the nucleus. Since the main interest 
of this paper is the extended emission, here we will 
focus our discussion on the nebular lines. 

For both 4C 41.17 and 4C 60.07 the oxygen line 
profiles closely resemble the bright inner parts of 
the Lya emission. Since self-absorption cannot be 
important in the strong forbidden lines, this sug- 
gests that the Lya line is giving useful informa- 
tion about the kinematic structure despite being 
subject to resonance scattering. This can only be 
possible if the interstellar gas is highly clumped. If 
this is true in the inner regions, it is probably also 
true in the outer regions. In any case, even if the 
line is being resonantly scattered in the outer re- 
gions, it would likely still provide a fair (although 
possibly velocity biased) measure of the neutral 
gas kinematics. 
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Fig. 4. — Similar to Fig. [3] but for the Lya emis- 
sion along the slits that lie more or less perpendic- 
ular to the radio axis of 4C 41.17. PA=19° (top), 
PA=22° (middle), and PA=171° (bottom). 

Figures El IU [6l and [lOl summarize radial ve- 
locities, FWHMs, and the run of relative sur- 
face brightness of the halos as a function of dis- 
tance from the radio core for different position 
angles. These figures show clear evidence for a 
distinction between disturbed and more quiescent 
regions. The inner regions have higher surface 
brightnesses, are characterized by large velocity 
dispersions (FWHM w 1500 kms" 1 ), and seem 
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to be embedded in low surface brightness regions 
with FWHMs of order w 500 kms" 1 . In g]we will 
discuss the general implications of these diagrams. 
First, we describe the individual sources. 

3.1. Notes on individual objects 

3.1.1. 4C 41.17: General Remarks 

iDev et all (1997) have discussed the brightest 
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RalaLrv's VelodV [kffis-') 
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T-Tt SE 5 r - 



part (2" x 1") of the extracted 1-D optical spec- 
trum of 4C 41.17 at PA = 76° in detail. They 
determined a redshift z = 3.79786 ±0.00024 based 
on the He II line. Furthermore they found evi- 
dence for stellar absorption lines, and low polar- 
ization indicating that a young stellar population 
contributes significantly to the rest-frame UV con- 
tinuum emission. 

The 2-D spectra along the radio axis (PA 76° 
plotted in fig. [3} show that the Lya, [O II], and 
[O III] emission line regions are very extended 
(over approximately 20", 10", and 6", respec- 
tively). The Lya and [O III] lines both show 
two separate components straddling the radio core 
with peak fluxes separated by about 3". Care- 
ful inspection reveals that these components are 
present also in the [O II] line, with the red part of 
the western component missing due to a skyline 
near A = 17880 A. The dip in between the emis- 
sion peaks is expected since the narrow-band Lya 
and radio image ov erlays showed the cor e to be 
highly obscured (see Reuland et al. ( 20031 )). 

The large extent of the nebula in the [O II], 
and [O III] emission lines shows immediately that 
the whole Lya, halo is enriched in heavy elements. 
This immediately discounts the possibility that 
the emission line halo represents in-falling pristine 
gas. This point will be discussed further below. 

A second item of note is that the velocity gra- 
dient in the halo and the velocity dispersions ob- 
served in the lines are comparable. This shows 
that the halo is a highly turbulent and disordered 
structure. 

The image overlays indicated that the position 
of the radio core is offset by 0.5" to the NE from 
the central dip in the emission. In the spectra, the 
position of the radio core is revealed by a narrow 
high velocity tail on the [O III] profile, extending 
out to -1500 kms -1 and clearly visible in figure 
[2J In the Lya profile, the nucleus is marked by a 
faint tail extending out to over +2000 kms^ 1 the 
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Fig. 5. — Similar to fig. [2] Grayscale represen- 
tations of the 2-D spectra of emission line halo 
around 4C 60.07. The top left and top right panel 
respectively show Lya and C IV emission along 
the filament with PA = 135°. The bottom left 
and right panel show Lya and [O II] A3727 along 
the radio axis with PA = 81°. 

negative velocity tail presumably being obscured 
by resonant scattering. 

3.1.2. 4C 41.17: Kinematics 



Dev et al.l (|1997l ) discerned the presence of both 
narrow and broad components in the central re- 
gions of the Lya and C III] emission lines of 
4C 41.17 with FWHMs of «500-600 kms" 1 and 
«1200-1300 kms -1 respectively. Figure [3] shows 
that these high velocity dispersions extend along 
the radio axis to the outer radio lobes. Beyond the 
limit of the radio lobes there is a break in velocity, 
the kinematics become more quiescent, and the 
velocity centroids change from being blue-shifted 
(—500 to —700 kms -1 along the Western lobe) to 
near systemic. Perpendicular to the radio axis (see 
Fig. [4j the kinematics follow a smooth gradient 
and the broad emission lines are found only in 
the central region of «3" (?»20kpc) wide, iden- 
tified by its sharp peak in peak surface brightness. 
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This bright emission is associated with the inner 
radio lobes, and has been successully ascribed to 
radiativ e shocks at the boun dary of the expanding 
cocoon (jBicknell et alj|2000h . 

Except for the high velocity tail originating 
close to the AGN and discussed above, the Oxy- 
gen lines show narrower line components than the 
Lya line, with FWHMs «600 kms -1 . The mea- 
sured velocities and velocity dispersions likely re- 
flect the true motions of the gas, whereas Lya 
is additionally broadened by resonant scattering 
into the damping wings. The measurements on 
the [O III] line are likely to be much more reliable 
than the [O II] line because the line is brighter, 
is less reddened, and is not so badly cut about 
by night sky lines. Furthermore, corrections to 
the FWHM do not need to be made (in con- 
tradistinction to the [O II] doublet), and finally, it 
can be observed closer to the nuclear region than 



, because of its higher critica l density (see 
di Serego Alighieri et al. 19971 ). 



[O II 
e-g-, 

Two very important results are found from the 
near-IR spectroscopy. First, as has been noted 
in the previous section, despite the obvious res- 
onant broadening, the velocity structure of Lya 
closely resembles that of [O III]. Secondly, there 
is [O II] emission in the velocity regime —200 to 
— 1000 kms -1 extending from the nucleus to 7-8" 
west of the nucleus, whereas the [O III] emission 
appears to be much more centrally concentrated. 
We will discuss this further in t]4.1l 

The Lya emission appears to show a strong ve- 
locity shear in its kinematically quiescent outer 
parts as yet undisturbed by the expansion of the 
radio source along the radio axis. By contrast, 
the kinematic profiles obtained with slit positions 
perpendicular to the radio axis (see fig. [4} show 
a fairly symmetric velocity distribution. This is 
indicative of an overall shear or rotation of the 
Rj 200 kpc diameter Lya halo about its major axis. 

The velocity structure of the bulk of the emis- 
sion line gas appears predominantly red-shifted in 
the Oxygen, H/3, and Lya lines. A similar red- 
shift is seen in the extend ed CO J=4— 3 emis- 
sion (|De Breuck et all [2005) , suggesting that the 
AGN may be offset from the systemic redshift of 
the galaxy. The central CO component is situ- 
ated at a relative velocity of —125 kms -1 . The 
location of this CO component coincides in spa- 
tial position and velocity with the Lya gap be- 



tween the "cloud" an d the galaxy as discussed in 
Reuland et al. ( 2003h . suggesting that the dense 
molecular gas is absorbing the Lya emission at 
this location. 

3.1.3. 4C 60.07: General Remarks 

The ESI spectrum of 4C 60.07 (PA = 135 °; 
fig. O shows a clear trace of continuum emission 
over approximately 3" which we associate with 
the position of the radio core and the unobscured 
optical nucleus. Comparison of this identification 
with the extent of the Lya filament extending to- 
wards to N-W narrow-band image (Fig. [1]) and 
which is visible in the lower part in the spectrum, 
shows that the nucleus has been located to within 
0.5 ". Based on the He II line at 7855.7 ± 1.1 A we 
infer a redshift z = 3.7887 ± 0.0007, in agreement 
with IRdttgering et al.1 (|1997h . 



3.1-4. 4C 60.01: Kinematics 

The 2-D Lya emission line (Fig. [6]) shows inter- 
esting structure: crescent shaped clouds surround 
a gap at the position of the radio core near the sys- 
temic velocity of the galaxy. The spatial structure 
was expected from the narrow-band image, but the 
kinematic substructure is unusual. A comparable 
lack of Lya emission near the systemic velocity of 
the galaxy has been found in a Lya galaxy associ- 
ated with the SSA-22 protocluster ( Wilman et al.l 
l2005h and the large Lya halo r ecently discovere d 
by the Spitzer Space Telescope ( Dev et al~ll2005h . 
It seems reasonable to attribute these velocity 
profiles to complex radiative transport effects of 
the resonantly scatter ed Lya emission in dense 
(dusty) media (e.g., iNeufeldl Il990l : lAhnl l2004h 
since these sources are strong submillimeter (rest- 
frame far-IR) emitter s ( Papadopoulos et"al] [2000: 
Chapman et al.ll200l[ ). However, the similarity of 
the profiles of the the Lya line and the C IV line 
in this region suggest, rather, an interaction of the 
radio lobe with a cocoon of material. 

The morphologies of the Lya and [O II] spec- 
tra at PA = 81 (fig [5]) are remarkably similar, 
providing further justification for the use of Lya 
as a useful tracer of the dynamics of these sys- 
tems, despite the problems of resonance scatter- 
ing Both lines have very high velocity FWHMs of 
~ 100 A or w 1600 — 1700 kms -1 , rather large for 
HzRGs. They are however smaller than for the 



8 



SE 5 



d 



-2000 



Relative Velocity (km s -1 ) 
-1000 1000 2000 



O 



-10 - 



NW 



T 



T 



T 



-5 - 




-L 



578 580 582 584 
Lambda (nm) 



586 



Fig. 6. — An expanded grayscale representation 
of the 2-D spectrum of 4C 60.07 centered at the 
redshifted Lya line showing the filament and the 
crescent shaped arcs surrounding a depression in 
the emission that is identified with the location 
of the both radio core and unobscured optical nu- 
cleus. This kinematical signature is probably the 
result of a shock interaction between the radio lobe 
and its cocoon. 

spectrum at PA = f35° which shows a FWHM 
of « 2600 kms -1 . All of this is suggestive of 
the strong radiative shocks produced in an interac- 
tion of a powerful radio source with a surrounding 
dense interstellar medium. 

The emission line profile across the Lya gap can 
be fit with two Gaussian components with central 
wavelenghts of 5839.9 ± 0.5 A and 5806.8 ± 0.4 A 
and FWHMs of 12.8 ± 0.5 A and 10.6 ± 0.3 A re- 
spectively. However, it is probably better to fit the 
whole bright region to a single expanding shell by 
a velocity ellipsoid, estimated by fitting an ellipse 
to the position : velocity curve obtained for the 
two peaks in line intensity at each spatial position. 
The fit obtained depends somewhat on both the 
position angle of the observation and on the line 
which is used to measure it, but all are consistent 



with an expansion velocity of 1100 ± 200 kms 1 . 

3.1.5. 4C 60.07: The Lya Filament 

Perhaps the most striking feature of the emis- 
sion line nebula surrounding 4C 60.07 is the ex- 
tended Lya filament. Imaging showed this to 
be of fairly constant surface brightness, sharply 
bounded on the NE side, and much more "fluffy" 
and rather ill-defined on the SW side close to the 
main emission region. The tip of the filament is 
cospatial with a small galaxy suggesting a causal 
connection. In Figure [8] we present the spectrum 
of this filament obtained with the ESI instrument. 

Figures [6] and [7] show that the filament is offset 
from the systemic velocity of the radio source by 
approximately —200 kms . The central wave- 
length fit by a single Gaussian yields A « 5818.1 ± 
0.3 A or —175 kms -1 relative to the systemic ve- 
locity at z = 3.789. The line has a FWHM of 
20 ± 1 A corresponding to a (deconvolved) velocity 
dispersion of « 300 kms -1 . The fact that there is 
no evidence for a strong (> 100 kms -1 ) velocity 
gradient across the filament seems to indicate that 
it must be moving in the plane of sky. 

Observations with mm-interferometers of the 
redshifted CO J=4-3 transition in 4C 60.07 have 
found two kinemati cally and spatially separat e 
gaseous reservoirs 



( Papadopoulos et al. 2000l ). 



Recently this has been confirmed by VLA obser- 
vatio ns of the CO 3—1— transition (jGreve et al 
120051 ) . An overlay of the L ya image with the mm - 
observations, see fig. 5 in iReuland et al.1 ( 2003 ). 



suggests that the filament extends in a direction 
perpendicular to the major axis of the molecular 
gas, reminiscent of the (polar ejection) morphol- 
ogy of a galactic superwind. The CO 3=4— 3 
observations showed a gas-rich narrow component 
with a FWHM of « 150 kms -1 at a velocity of 
—224 kms -1 to the galaxy. This same compo- 
nent has now been independently dectected by 
its CO 3=1— emission at a relative velocity of 
-220 ± 40 kms" 1 with a FWHM of 165 kms" 1 . 

As is the case for 4C 41.17, the optical spectro- 
scopic observations show similar velocities for the 
molecular and emission line gas. This corroborates 
the idea of a superwind, as it provides a dynam- 
ical connection between the filament and the gas 
fueling the nuclear AGN. If this is the case, the 
source of the ionization may be either shocks or 
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Fig. 7.— Similar to fig. [3] but for 4C 60.07. Lya 
(top) and [O II] (middle) emission along PA=81 
and Lya emission along PA= 1 35 (bottom). 

photoionization from the central nucleus. 

The superwind interpretation of the Ly— a fila- 
ment does however present certain difficulties. In 
particular, it is hard to understand how this fea- 
ture, over 70 kpc in length, manages to remain so 
well-collimated and kinematically quiescent over 
such a distance if it were related to a high- speed 
outflow (albeit in the plane of the sky). However, 
the filter cut-off may play an important role here, 
since the the observed velocity of the filament is 



very close to the filter cut off. Thus, it may give 
the appearance of a sharp boundary, but may not 
be one in actuality. Possibly a slit placed over 
"empty sky" next to the filament would reveal 
other portions of the filament at different veloci- 
ties. As alternative possibility may be that it rep- 
resents an accretion filament or a tidally stripped 
gas stream extending from the active galaxy. In 
this case, the both the quiescence of its velocity 
field, and the similarity of the systemic veloci- 
ties of the filament and of the molecular reservoir 
would both find a more natural explanation. Pho- 
toionization by the central AGN would then be the 
most probable cause of the ionization of the fila- 
ment. On balance, this accretion or tidal stripping 
scenario seems to provide the most convincing ex- 
planation for the filament. 

The galaxy at the tip of the filament is most 
likely to be a chance superposition. Figure [8] shows 
the 2-D and 1-D spectra of this galaxy (Compan- 
ion # 1), and a second companion galaxy (Com- 
panion #2) which also fell on the slit. The spec- 
trum of Companion #1 is very similar to, but 
fainter than Companion =^2. The brighter galaxy 
shows an emission feature, identified as the [O II] 
line and a break near 7000 A, while the overall 
shape of the spectrum for the fainter galaxy is sim- 
ilar but without an evident break. Identifying the 
line at 7047.0A with redshifted [O II] A3727 yields 
a redshift of z — 0.891 for the galaxy. This implies 
that it is a foreground object and not causally con- 
nected to the filament, despite its suggestive loca- 
tion. 

3.1.6. B2 0902+34: General Remarks 

The optical spectrum obtained for B2 0902+34 
is more sensitive than previous observations re- 
ported in the literature (Lilly 1988; Martin- 
Mirones et al. 1995). We have detected Lya 
emission over an extent of ~ 10" with a complex 
multi-component spatial and velocity structure 
(cf. fig [9]) and have detected extended (3") con- 
tinuum emission from the core over the entire 
wavelength range with an almost constant surface 
brightness. The Lya profile shows galaxy- wide 
blue-shifted absorption by neutral hydrogen. 

As is evident from Tableland Figure [H the in- 
ferred redshifts depend on both the emission line 
and the position of the aperture used. Much of 
this is due to the absorption optical depth effects 
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tion of the galaxy might suggest a physical associ- 
ation with the Lya filament of 4C 60.07. However, 
the spectrum of companion galaxy #1 bears close 
resemblance to the brighter foreground galaxy #2 
shown in panels (b) and (c), respectively. The 
circles in panel (a) indicate emission lines which 
are probably redshifted [O II] A3727 in both fore- 
ground galaxies, identified in the spectra of panels 
(b) and (c). 

in the line. The He II line at 7199.35 ±0.5 A, ex- 
tracted in a 2" wide aperture centered at the radio 
core, is probably most representative of the sys- 
temic velocity of the galaxy. From this we infer 
a redshift z = 3.3886 ± 0.0003. Th is is slight!; 
less than that reported previousl y by iLillv ( 198 
and iMartin-Mirones et al. ( 19951) . They inferred 
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Fig. 9.— Similar to fig. H but for B2 0902+34. 
The top left panel shows the line emission profile 
for Lya, top right C IV, bottom left He II, and 
bottom right [O II]. 

z = 3.395 and z = 3.391 respectively. How- 
ever, these reported redshifts are entirely consis- 
tent with redshift estimates we derive using other 
lines in the same aperture, emphasizing the diffi- 
culty of assigning single redshifts to such complex 
and self-absorbed systems. 

3.1.7. B2 0902+34: Kinematics 

The He II line is only a few hundred kms 
wide and may be slightly extended, with its po- 
sition of maximum surface brightness centered at 
the position of the radio nucleus. We can therefore 
argue that it arises from photo-ionization of the 
dense gas surrounding the nucleus. The [O II] line 
is broad and spatially extended, but is too faint to 
allow a detailed comparison with the kinematics 
of the other lines. 

Given that B2 0902+34 shows associated H I 
absorption against the radio continuum it is im- 
portant to see whether this is present also in 
the Lya profile. This emission line shows a self- 
absorption feature on the blue side of the profile 
centered at A = 5331.58 (z — 3.3857), and which 
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can be fitted with a Voigt profile with a column 
density of 1.8 x 10 14 cm -2 and Doppler param- 
eter b = 195 ± 11 kms -1 . The inferred blue 
shift is « 100 kms" 1 . This cannot be the same 
gas that causes H I absorption against the radio 
contin uum, since iBriggs. Sorar, &: Taramopoulosl 
(|l993l ) and |Cody_jLBraun| (|2003j) infer a redshift 
z abs = 3.3962 with N H = 3 x 10 21 cm" 2 and 
FWHM = 120 krns" 1 for the absorbing gas. It is 
probable that the dense gas giving the H I absorp- 
tion lies close to the nucleus, and any absorption 
that this may have produced in the Lya line is 
completely veiled by the frequency re-distribution 
produced by multiple scattering events. 

The Lya absorption can possibly be identified 
with a galaxy-wide outflow of material. In this re- 
spect, it is v ery similar to the gal axy- wide outflow 
identified bv lWilman et alj (|2005l ) in a Lya galaxy 
associated with the SSA-22 protocluster. We will 
return to this point in the discussion. 

The underlying Lya emission profile revealed 
by fig. [9] is very interesting, as it shows evidence 
of a strong velocity shear amounting to nearly 
1000 kins" 1 across the central arc sec. The up- 
per (NE) component has the typical triangular 
shape resulting from spatially extended absorption 
that is blue-shift ed relative to the systemic veloc- 
ity (cf . iDev 19991 ). and it is terminated by the self- 
absorption on the blue side. The lower (SW) com- 
ponent may well have an underlying shape which 
is the symmetrically reversed version of the NE 
component, but the self- absorption makes this dif- 
ficult to trace, and the line reaches its maximum 
brightness at a relative velocity of —500 kms -1 , as 
it emerges from under the absorption. Because of 
the self absorption, a gaussian fit to the line profile 
is not valid, and yields spuriously large values, as 
can be seen in fig. [TO] 

The C IV A1549 emission is both spatially ex- 
tended and broad (« 1000 km s -1 ) and is centered 
on the nucleus. The line broadening in the nu- 
cleus is probably, once again, the result of a direct 
interaction between the relativistic outflo w asso- 
ciate d with the peculiar radio source (see Carillil 
19951 . for details) and its dense cocoon of interstel- 
lar medium. 



4. Discussion 

4.1. Emission Line Diagnostics 

The relative intensities of emission lines are, in 
principle, powerful diagnostic tools for studying 
the excitation mechanisms, metallicities and phys- 
ical conditions in the emission line gas associated 
with AGN. Diagnostic diagrams for optical line 
ratios have been most extensively used to probe 
the gas in nearby active galaxies, where evidence 
for both jet- and accretion- powered shocks and 
for photoionizat ion by the central AGN ha s been 
adduced (e.g.. iBicknell. Dopita & O'Deal Il997t 



Dopita et al.ll997t iGroves. Dopita &: Sutherland 
l2004al lbT). Consequently these relationships have 
been best calibrated for rest-frame optical lines. 
However, line diagnostic diagrams applicable to 
both shock- and photo-ioni zed AGN have been de- 
veloped for use in the UV (lAllen. Dopita k. Tsvetanovl 



pr 

19981 : IGroves. Dopita & Sutherland! l2004bD . 



Until recently, for most HzRGs (|De Breuck et al 



2000) or for the 3C radio galaxies (Best, Rottgering, 



& Longair 2000a,b; Inskip et al. 2002a,b) only the 
rest-frame UV lines have been available. The opti- 
cal line diagnostics have previously not been com- 
monly used because of calibration difficulties and 
the limited sensitivity of near-IR spectrographs on 
telescopes of the 2-4 m class. 

We would like to be able to use both UV and 
optical emission line diagnostics to determine the 
mechanism that is responsible for ionizing these 
extended emission line nebulae. The most likely 
candidates are: photoionization by radiation from 
an AGN or stars, shock heating, and shocks with 
precursors that photoionize the region ahead of 
the shock by radiation from the gas. Studies 
of z ~ 1 radio galaxies have shown that the 
dominant ionization mechanism may depend on 
the evolutionary state of the rad io source (e.g., 
Best, Rottgering &: Longair 200 0a, b: In skip et al 



2OO2a0). Generally, it was found that the emis 



sion line gas of small (i.e. young) sources is shock 
ionized, but as the radio source expands beyond 
the host galaxy, interactions with the gas decrease 
and photoionization by the AGN takes over. As 



De Breuck et al.1 (|2000h noted, for HzRGs the sit- 



uation may be more complex. 

Before we investigate these line diagnostics, we 
shall first use global energy requirements to con- 
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strain important source parameters on the basis 
of both photoionization and shock models. The 
total Lya fluxes of the extended emission line 
regions are of order 10~ 14 ergs -1 cm" 2 , yielding 



lumin osities LLya « lO^ergs" 1 (jReuland et al 
(l2003h . If the source of ionization is an active 



nucleus, then the luminosity of the source derives 
either from the UV photons which it produces, or 
from the mechanical energy of the jets which it 
powers. These luminosities can be derived from 
the theoretical ratio of the Ly— a flux to the to- 
tal flux in both shock and photoionization mod- 
els. To do this we use th e shock models from 
Dopita fc Sutherland! ( 1996h . and assume that the 
PdV work inmplied by this luminosity is the me- 
chanical energy flux of t he jet. For the justification 
of this assumption, see iBicknell. Dopita fc O'Dea 
(|1997l ) . For photoionization, we use the radiation- 
pressure dominated photoionization models from 
Groves, Dopita, & Sutherland (2004a,b), assum- 
ing a power-law spectrum with index a — —1.4 
(/„ w v a ). From these we derive either an ioniz- 
ing radiative or a mechanical energy flux of at least 
1.1 x 10 46 ergs _1 . This flux is consistent with a 
typical luminous embedded QSO. 

Alternatively, some of the ionization may be 
due to young stars such as are found in the nu- 
merous clumpy co mponents of the Southe rn re- 
gion of 4C 41.17 (jvan Breugel et alJ fl999h . and 
whose formation has probably been triggered in 
the dense cocoon surrounding the cent ral radio 
jet (|Dev et all 1 199^ IBicknell et al.ll2000h . In this 
case, assuming a stellar radiation field with an 
effectiv e temperature of 42,000 K, we infer (us- 
ing the Kennicuttl ( 19981 ) calibration) that a star 
formation rate of at least 3000 M Q yr -1 would 
be needed to ionize the halo. This star forma- 
tion rate is remarkably close to w hat is inferred 



from rest-frame far-IR ob servations ( Du nlop et al 



1994t IStevens et al.ll2003h . However, this does not 
necessarily imply that stellar photoionization is 
the dominant ionization mechanism. Photoion- 
ization by stars usually results in Lya equiv- 
alent widths less than 240 A (a ssuming stellar 
populations at solar metallicity; ICharlot fc Fall 
19931 ) and would seem inconsistent wit h the large 



equiv alent widths observered (« 500 A; iDev et al 
1997). The strength of the C IV lines and other 
high-excitation species also militates against a 
purely stellar origin of the ionizing flux. Indeed, 



Bic knell et al. I (|2000h have used the emission spec- 
trum to infer that shocks are the dominant source 
of ionization in this region, and have used the 
flux to determine the PdV work being done by 
the relativistic jet. However, away from the ra- 
dio axis photoionization by stars may well be- 
come more important. 4C 41.17 is a young sys- 
tem and recent modeling suggests that any pri- 
mordial stellar populations could be very efficient 
in emitting significant amounts of ionizing flux 
while re maining virtually undetectable in the opt i- 
cal (c.f., iFosburv et al.ll2003tlPa~nagia et al.ll2003h . 

Using the observed line ratios, better con- 
straints on the ionizing process can be derived. 
Dev et al. I (|l997l) determined the following ratios 



for the narrow rest-frame UV lines of 4C 41.17: 
C III]/C IV w 0.14 (w 0.7 for total fluxes) and 
C IV/He II w 2.4. The latter is close to the maxi- 
mum of 3.1 predicted by photoionization models. 
Although the narrow emission line strengths of 
C III], C IV, and He II can be matched by a 
simple nuclear AGN photoionization model with 
solar metallicity clouds, a high ionization param- 
eter U < 0.1, and a po wer-law ion i zing s ource 



with index a = -1.5, IDev et all (|l997l) and 



Bicknell et al. ( 2000h argue that the dominating 
ionizing mechanism is shocks rather than pho- 
toionization the nucleus. 

Figure [Tl] shows that a combination of rest- 
frame UV (C IV A1549/He II A1640) and rest- 
frame optical ([O III] A5007/H/5) emission line 
flux ratios can help to separate pure photoion- 
ization from shock dominated mechanisms. This 
diagram was derived assuming solar metallicity 
gas and is relatively insensitive to the effects of 
dust extinction as the ratio s are dete r mined from 
lines close in wavelength. Ilwamur 3 (|2003h con- 
ducted a rest-frame UV-optical emission line study 
of 15 radio galaxies with 2 < z < 2.6. They 
found that there is a range in observed line-ratios 
suggesting that some objects are best explained 
with photoionization of low metallicity gas while 
others are consistent with the s hock+precursor 



model . Carson et al. < 200ll ) and Maxficld et al 



(J2002J) found evidence for changing line ratios 
within sources, suggesting that the dominant ion- 
ization process are depends on the region of inter- 
est. 

Based on previously published values of [O III]/ 
H/3 w 3.4 and 2.8 for 4C 41.17 and B2 0902+34 re- 
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Fig. 10.— Similar to fig. [3] but for B2 0902+34. 
Lya (top) and [O II] (bottom) both along 
PA=69°. 



spectively (jEales &: Rawling slll993h it seemed that 
pure photoionization models could be ruled out. 
However the H/3 detections were marginal, and us- 
ing the total integrated line fluxes for 4C 41.17 
from the present study we derive ratios of [O III] 
A5007/H/3 «9.6 for the central region and [O III] 
A5007/H/3 w 11.8 for the spatially integrated spec- 
trum that are much higher and are consistent with 
both pure photoionization and shock+precursor 
models. The point derived for the spatially in- 
tegrated spectrum of 4C 41.17 is also shown on 
fig. [HJ This appears to show that the excitation 
of 4C 41.17 is similar to other radio galaxies. It 
cannot be due predominantly to stars, for which 
the [O III] /H/3 would be of order three and for 
which both C IV A 1549 and He II would be very 
faint. 

The extended regions of 4C 41.17 have a much 
lower excitation, however. As discussed in <j3 .1.11 
for 4C 41.17, [O II] emission was detected as far 
as w 60kpc from the nucleus. Figure [T2l shows the 
velocity integrated relative intensities of Lya and 
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Fig. 11. — A UV-o ptical line diagnostic dia- 
gram adapted from en. Dopita fc Tsvetanovi 
(|l998h . The ratios [O III] /H/3 and C IV/He II 
are used to separate shock (lower grid), 
shock+precursor (right hand side) and pho- 
toionization (top left) models for solar metallicity 
gas. The shock models are marked according 
to the shock velocity and magnetic parameter, 
and the photoionization models according to 
the ionization parameter, U, and the power-law 
index of the photoionizing radiation field (a). 
The points represent nearby z » 1 Galaxies, and 
the accretion disk of M87 in explicitly identified. 
The large point with error bars is for the nucleus 
of 4C 41.17. This proves that the excitation of 
4C 41.17 is similar to the other radio galaxies, 
and like them is shock-dominated. 



[O II] as a function of distance from the nucleus for 
a red-shifted (+200 to +700 kms -1 ) and a blue- 
shifted (-200 to -1000 kms -1 ) velocity interval. 
The red-shifted [O II] flux closely follows that of 
the Lya , while the blue-shifted [O II] shows a 
relative enhancement over the range 3-8" W of 
the nucleus. In this region no evidence for [O III] 
is found (Fig. [2]), whereas the [O III] emission 
in the nuclear region is much brighter than the 
[O II]. This shows that [O III]/[0 II] ratio is much 
lower («l-2) in the outer regions than in the cen- 
ter («3-4) implying that a different mechanism is 
responsible for the extended [O II] emission. The 
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Fig. 12. — Top: Normalized surface brightness 
profiles of the [O II] (solid line) and Lya (dashed 
line) emission along the inner radio axis and south- 
west filament of 4C 41.17. These are integrated in 
the velocity range +200 to +700 km s _1 and clearly 
show the spatially-extended [O II] emission. The 
spatial zero-point corresponds to the position of 
the radio core. Bottom: Similar to top panel but 
for the velocity range —200 to — lOOOkms -1 . The 
[O II] emission is detected out to w 8" (~ 60 kpc) 
west of the nucleus, where it is blue-shifted by 
w — 600 km s -1 relative to the He II line. The pro- 
jected distance of the south-west and north-east 
radio lobes along the slit direction are indicated 
(dotted lines). The +200 to -200kms" 1 range 
for [O II] is confused by near-infrared sky lines 
and is therefore not shown. 



[O III]/[0 II] ratio is known to be sensitive to the 
ionization parameter in the nebula. Because they 
have higher ionization parameters and also have 
central stars with higher effective temperatures for 
a given age, the low-metallicity H II regions are 
characterized by higher [O III]/[0 II] ratios. Since 
this ratio falls to very low values in the outer parts 
of 4C 41.17 this proves that Population III stars 
with their very high effective temperatures cannot 
be responsible for the excitation of the outer re- 
gions. Excitation by shocks is a more likely ioniz- 
ing source, since the emission is extended along the 
radio-axis and both the Lya and [O II] are blue- 
shifted by w 600 kms -1 . Furthermore, the peak 
in the [O II] emission on the W side of 4C 41.17 
corresponds to the point at which the velocity dis- 
persion drops steeply and boundary of the outer 



radio lobe. This suggests that the [O II] is ex- 
cited by shocks in an outward-propagating cocoon 
around the radio lobe. We note here that [O II] 
emission is particularly stro ng in shocks. This cor- 
roborates the suggest ions bv lDev et al.1 ( 1997t ) and 
Bick nell et al. | (|2000l ) that the Lya emission, at 



least in this direction, is not due to scattering but 
that the nebula must be locally ionized by s hocks 



related to radio source. IScharf et alJ ([2003) 



pro- 



posed that the escape of X-ray emission from the 
shocked region could be responsible for ionizing 
the extended halo beyond the the radio lobes. 
Fro m the jet shock spectrum, iBicknell et al 



(2000) inferred that the metallicity of the star- 
forming cocoon around the radio jet was about 
the same as for the LMC; a half to a third of so- 
lar. For the nucleus itself, N V is often used as 
a tracer of the chemical abundance, because it is 
rather insensitive to shock ionization, density fluc- 
tuations, ionization parameter, and as a secondary 
nucleosynthesis element has a quadratic depen- 
dence on the metallicity jVillar-Martm et aljl999l 



iGroves. Dopita fc Sutherland! l2004ah . This would 
imply that that would only be detected in highly 
chemically enriched sourc es, and not significant 
in the majority of HzR G s dRottgering et al1ll997t 
De Breuck etaH l200d) . IVernet et all (|2001l ) pre- 



sented a metallicity sequence for HzRGs that is 
simil ar to the one for quasars ( Hamann fc Ferlandl 
1993). 4C 41.17 was included in their study and a 
metallicity of Z ks 1.3 Z Q was i nferred for the ce n- 
tral region (aperture of 2" x l" |Pev et ailll997l) . 

Thus, the nuclear and jet cocoon spectra of 
4C 41.17 argue for a fair degree of chemical enrich- 
ment, perhaps more in the nucleus than in the jet 
region. However, the halo region is by no means 
primordial in its abundances. Although we can- 
not directly infer the actual metallicity, we can say 
that the detection of strong [O II] emission out to 
w 60 kpc from the nucleus demonstrates that the 
halo has been chemically enriched throughout its 
volume, at least in the illuminated regions. 

A similar conclusion is also inferred for 4C 60.07. 
The very similar velocity morphologies and spa- 
tial extent of the Lya and [O II] spectra at PA 
= 81° seen in figure [5] once again shows that the 
enriched gas is distributed throughout the Lya 
halo. The strongly- C IV emitting shell in the 
nuclear vicinity shows that this galaxy too is ex- 
cited primarily by the AGN rather than by star 
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formation. The expanding shell morphology of 
both Lya and C IV further argues that 4C 60.07, 
like 4C 41.17, is excited mostly by shocks in the 
dense expanding cocoon of interstellar gas pushed 
by the relativistic jets from the nucleus. 

4.2. Galaxy Mass estimates 

The HzRGs are known to represent some of 
the most massive galaxies galaxies formed in the 
early universe. An estimate of their mass, by any 
means, is important to constrain the A-CDM mod- 
els of galaxy formation, since the formation of such 
massive galaxies so early on might present some- 
thing of an issue for the theory to resolve. We 
therefore provide these mass estimates in this sec- 
tion. 

4- 2.1. Luminosity based masses 

van Breugel et al. (1999) estimated the mass of 



4C 41.17 using the integrated rest-frame UV light. 
They found it to be a ver y massive galaxy with 
L uv « 17 L*. Similarly, I Graham et "all (| 19941 ) 
used near-IR imaging with the line-free K$ fil- 
ter to estimate the mass of the galaxy directly. 
In a 4" diameter aperture they found K$ = 
19.6 ± 0.6 mag. This converts in their cosmol- 
ogy to an absolute magnitude of 4C 41.17 to be 
M B = -23.0, or about 20 L*. 

A lower mass limit for the ionized gas in the 
halo can be obtained from photoionization mod- 
elling of the total Lya flux, neglecting the ef- 
fect of dust absorption, and assuming that the 
ISM is fully ionized. From this we derive an 
average thermal electron density in the halo of 

_JL 

n e ,haio = 0.1 / v 2 cm -3 where / v is the volume 

filling factor of the Lya emitting gas. This yields 

i 

a mass of gas of M ha i ~ 6 x 10 12 / v 2 M . For a 
volume filling factor / v ~ 1 the inferred amount 
of gas is substantial and comparable to the broad- 
band mass estimates of the host galaxy and that of 
fully formed, massive elliptical galaxies in the lo- 
cal Universe. However, the fact that the Lya line 
appears to provide a good dynamical mass tracer, 
and is not grossly affected by self-absorption ar- 
gues in favor of f v << 1, and consequently, for 
a much lower total mass. The molecular com- 
ponent of 4C 41 . 17 ha s recently been traced by 
|Pe Breuck et al.l ( 2005 ). who have detected the 
CO(4-3) transition using the IRAM interferome- 



ter. They find two massive CO clouds « 
6 x 10 10 M©) which coincide with two different 
dark lanes located close to the nucleus in the deep 
Lya image. This is similar to the CO content 
of the sub-millime tre galaxies, as measured by 



Greve et al.l (|2005l ) (M CO « 3 x 10 10 M ), sug- 



gesting once again that the HzRGs and the sub- 
millimeter galaxies are drawn from similar popu- 
lations of galaxies, differing only in the activity of 
their AGN. 



2.2. Dynamical mass estimates 
van Oiik et all (|l997h and IVillar-Martin et al 



(|2003l ) have presented mass estimates for the more 
quiescent parts of the emission line regions in their 
studies. There are two methods of estimating the 
dynamical masses. First, one can assume that the 
halos consist of gas that has settled in rotating 
disks. Secondly, the halos can be envisaged as 
consisting of virialized clumps, which have velocity 
dispersions that balance the gravitational forces. 

In the case of rotating disks, the mass can be es- 
timated by measuring the velocity shear across the 
halo and using: M|°^ = RV 2 /Gsm 2 i with R the 
radius of the disk, V half the amplitude of the ro- 
tation curve, and i the inclinati on of the disk with 
respec t to the plane of sky. IVillar -Martin et al 



(|20031 ) found evidence for rotation in about s» 50% 
of the objects they looked at and inferred masses 
of order M™* n x sin 2 i w 0.3 - 3 x 10 12 M . 

Looking at our measured relative velocity dis- 
tributions for 4C 41.17, we infer a velocity shear 
of « 300 kms" 1 (+100 kms -1 to -200 kms^ 1 ) at 
PA = 76 °, over a distance of R = 20 "/2 = 77kpc. 
From this we infer a dynamical mass of about 
M\f n x sin 2 i w 4 x 10 11 M©. The velocity pro- 
files at position angles PA =19 °, and 21 ° show no 
signs of a shear, consistent with their approximate 
alignment along the axis of rotation. 

The observations for 4C 60.07 does not show 
an indication of being settled in a disk. Rather, 
it could be supported against gravitation by ve- 
locity dispersions of the cloudlets in the halo. 
For this scenario the dynamical mass is given by: 
M™ n = 5RVg/G, with V R the radial velocity dis- 
persion of the clouds. Using this relation we de- 
rive dynamical masses for the outer regions of 
M dyn = 5 - 10 x 10 11 M Q . For the inner regions, 
where the signal-to-noise ratio is better we find 
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lower mass of M™ n w 2 x 10 11 M Q . 

The Lya emission profile of B2 0902+34 re- 
vealed by fig. [9] suggests a strong nuclear ve- 
locity shear. If this is interpreted as a rota- 
tion, then it would imply a rotational velocity of 
w 550 kins" 1 over a distance of 3" or a rotation 
velocity of rs 275 kms -1 at i2 « 12 kpc. This 
implies a dynamical mass within this core region 
of M|°^ x sin 2 i « 2 X 10 u M Q . However, as dis- 
cussed in section 13". 1.71 this velocity gradient could 
equally well be due to the expansion of a shell of 
interstellar gas driven by the a directed bipolar 
outflow from the AGN with a velocity w CX p- 

5. Outflows 

Absorbing gas associated with lower velocity 
outflows has now been observed in several HzRGs, 
almost all of which show asymmetric Lya pro- 
hies suggesting the presence of blue-shifted ab- 
sorbing gas which appears to be spatially ex- 
tended over the entire emission line region (s ee 





van Oiik ct al 


1997 


2003; 


Wilman et al.l 


2oo4. 



for outflows at high redshifts exists also for other 



presumably less massive , galaxies (jPettini et al 



20011: iDawson e t al 2002) 



There is an ongoing debate whether massive 
(> L*) galaxie s or less massiv e (< L*) galaxies en- 
rich t he ICM. iMartinl (|l999l ) and iHeckman eTal 



(2000) found that outflow speeds, y exp , are largely 
independent of galaxy mass. This would imply 
that smaller galaxies are more efficient at eject- 
ing enriched ma terial to large radii. Specifically, 
Heckmanl ( 2002 ') claim that for massive galaxies 
the metals will not escape the deep potential wells. 
However, recent modelling suggests that the bulk 
of the metals in cluster s are produced by a nd 
brighter galaxies (e.g., Nagashima et al.l 2004 ). 

For all the three HzRGs we have observed, we 
found features that can be interpreted as symp- 
tomatic of outflows. The extended optical fila- 
ment of 4C 60.07 is consistent with an outflow 
in the plane of sky. 4C 60.07 also shows cres- 
cent shaped arcs around the radio core and un- 
obscured optical nucleus that are consistent with 
an expanding cocoon; see fig. [Sj The blue-shifted 
carbon-rich component of B2 0902+34 seems to 
point towards an outflow of enriched material with 
an outflow velocity, v cxp of up to 1000 kms . 



The near-infrared [O II] spectroscopy (Figs. [U \3\ 
and [T2"|) and the earl ier optical Lya s pectroscopy 
of 4C 41.17 (Fig. El iDev et al.lfl997h . show that 
both emission lines exhibit large blue-shifted ra- 
dial velocities rj 600 — 900 kms" 1 seen in pro- 
jection along the radio axis. In particular, the 
gas is very disturbed along the south-west fila- 
ment, with Lya velocity widths ranging up to 
Av F ~ 900 - 1600 kms" 1 . Beyond the radio 
hotspot, the velocity and velocity widths decrease 
abruptly. The kinematics, its radial filamentary 
structure, and chemical enrichment of the gas all 
indicate a process of entrainment of material away 
from the central regions by the radio jet. In 
this scenario, the optical filament represents the 
shocked radiative cocoon of the radio source. 

For a galaxy mass of M ga i ~ lx 10 12 M Q within 
a radius of R ~ 40kpc, figures that seem to 
be typical for these HzRGs, the escape velocity 
v csc = (2GA//i? ga i) 1 / 2 is of order - 450 km s" 1 . 
The velocity gradients found from the spectro- 
scopic observations range a factor 1.5-2.5 larger 
than this, and implies that the ejection of the ISM 
originating from close to the nucleus of such galax- 
ies is a viable scheme to both chemically enrich the 
inter-galactic medium and to terminate the epoch 
of star formation in these massive galaxies. 

While it may still be true that purely starburst 
powered super galactic winds may not be suffi- 
ciently energetic for enriched nuclear material to 
escape from the galaxian potentials, the additional 
driving force produced by the energy content of 
the relativistic lobes produced by the central AGN 
seem quite able to blow matter out of the poten- 
tial of the galaxy. This is direct evidence of the 
op eration of the type o f feedback process proposed 
bv lSilk fc Reesl (| 19981 ). In this, the growth of the 
super-massive central black hole is terminated by 
the ejection of the reservoir of interstellar gas, once 
the mass of the black hole exceeds a critical value. 
In their model, this critical value is determined 
by the relationship between the radiation pres- 
sure and the depth of the stellar potential, and it 
is able to provide a qu alitative description o f the 
Magorrian relationship (jMagorrian et al . 1998). or 
the relationship between the mass of the central 
black hole and the velocity dispersion of the bulge 
( Magorrian et alJll998t iFerrarese fc MerrittJl200Cl: 
Gebharddl2000l) . 



The observations presented here suggest that 
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this model needs to be modified somewhat. 
Rather than being driven by radiation pressure, 
the dynamical evidence is that they are driven by 
the PdV work done on the galactic medium by 
the relativistic jets, and that the ejection is occur- 
ring in the faster moving material in the shocked 
cocoon around the radio jet. Such a scheme is 
also su ggested by the hydrodynamical simula - 
tions of ISpringel. Pi Matteo. fc Hernquistl (|2005l ). 
However, rather than being ejected isotropically 
in a disk-like galaxy, as in their model, in the 
HzRGs directional jets are being driven into a 
more spheroidal distribution of galactic interstel- 
lar gas. Given that 4C 41.17 has clear evidence 
of large amounts of shock-triggered star formation 
in this shocked cocoon, we can conclude that, in 
the HzRGs we are witnessing the fireworks that 
terminate the epoch of star formation in these 
massive galaxies - the last hurrah before these 
galaxies become massive "red and dead" Ellipti- 
cal galaxies that we find today in the cores of the 
most massive clusters. 

The multiple component, asymmetric, and 
twisted radio structure of 4C 41.17 suggests that 
during this phase, the central SMBH has ex- 
perienced multiple peri ods of radio sour ce ac- 
tivity and precession. (Ivan Breugel et ab 



1999; 



Steinbring. Crampton &: Hutchingsl 120021 ). The 
precessing radio source will gradually evacuate 
the central region of the galaxy of its interstel- 
lar medium. The time scales for star formation 
and radi o source activity in 4C 41.17 a r e very 
similar dChambers. Milev fc van BreugeJ 4990; 
Bicknell et al. 120001 ) and comparable to that for 
transporting the [O II] gas along the filament out 
to the vicinity of the south-west hotspot (« 7x 10 7 
yrs at w 900 kms -1 ). It suggests an overall pic- 
ture where the initial growth of the galaxy through 
merging, the feeding and growth of the central 
black hole, the growth of the stellar component, 
and finally triggering of SMBH the outflow, and 
its associated starburst activity are all closely cou- 
pled processes. 

6. Conclusions 

Extended Lya halos are observed both in 
HzRGs and in the sub- millimetre galaxies. In- 
deed, these objects may be drawn from the same 
underlying population, with the distinction that 



the HzRGs represent the radio-loud and active 
phase of the galaxian evolution. This seems sup- 
ported by the observation that the largest radio- 
quiet Lya halos have similar sizes and only slightly 
lower luminosities as the HzRGs but lack their 
large multi-component continuum structures. 

Our long slit optical spectra of 4C 41.17 , 
4C 60.07 and B2 0902+34 have shown that their 
Lya halos exhibit disturbed kinematics, with 
broad lines, large velocity shears, and, in some 
cases, expanding shells associated with the radio 
lobes. This clearly demonstrates that the relativis- 
tic jets are driving strong shocks into the galactic 
medium. 

In 4C 41.17 the near-IR spectra reveal very ex- 
tended [O II] and [O III] emission distributed along 
the radio source axis, as far as 60 kpc from the 
nucleus. This provides direct proof that the Lya 
halos are both chemically enriched by star forma- 
tion and ionized throughout the majority of their 
volume. The hypothesis that the Lya is due to 
scattering off H I clouds is disproved. Likewise, 
the idea that the halo gas is primarily composed 
of chemically pristine "primordial" gas is also dis- 
proved. However, we cannot be certain that the 
gas in directions perpendicular to the radio axis 
contains a large primordial component, because 
the sensitivity of our spectroscopic observations 
is insufficient to detect Oxygen lines at the low 
surface brightness expected. To settle this point 
would require very deep near-IR spectroscopy of 
the halo gas outside the radio sources, well be- 
yond the hotspots and orthogonal to the radio 
axes. This would be a very difficult observation 
to make. 

Our observations of HzRGs also help cast ob- 
servational light on the origin of the correlation 
found in galaxies between the ste llar velocity dis- 
persion and the black hole mass (lGebharddl200rl 
Ferrarese fc Merrittl l2000h . We find clear evi- 



dence of outflows of chemically enriched gas as- 
sociated with the jets, and velocity dispersions 
in the expanding cocoons around the jets which 
probably exceed the escape velocity. In addi- 
tion, B2 0902+34 shows Lya absorption can be 
identified with a galaxy-wide outflow of material, 
similar to the galax y-wide outflow identified by 
Wilman et aD (|2005h . This demonstrates that tin 



the HzRGs the central black hole has grown suffi- 
ciently to be able to profoundly modify its sur- 
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rounding galactic medium, pushing it into out- 
flow, and (in the case of 4C 41.17) triggering large 
quantities of star formation in the shocked galactic 
medium cocooning the radio lobes. Furthermore, 
the precession of the radio lobes seen in 4C 41.17 
will ensure that eventually, all of the interstellar 
medium will be shocked and either be converted 
into stars, or else ejected from the galaxy. 

These observations ca n be understood in terms 
of the general picture of Silk &: Reesl (1998) that 
feedback of the black hole on its host galaxy even- 
tually limit s the growth of the black hole. However 
unlike the Silk k, Reesl ( 19981 ) concept, the feed- 
back process is not primarily radiation pressure 
but is the mechanical energy input delivered by 
the relativistic jets. Thus, in the HzRGs we may 
be observing the moment where galaxy collapse 
gives way to mass ejection, a newly born galaxy is 
revealed. We may speculate that this is the defin- 
ing moment where the such galaxies enjoy one last 
violent burst of shock-induced star formation be- 
fore beginning their evolution to become the "red 
and dead" massive Elliptical galaxies we see in our 
local universe. 
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Table 1 
Radio galaxy sample 



Source 


KA (J 2000) 


DEC (J 2000) 


z 


B2 0902+34 


09 05 30.11 


34 07 55.9 


3.389 


4C 60.07 


05 12 55.17 


60 30 51.1 


3.789 


4C 41.17 


06 50 52.14 


41 30 30.7 


3.798 



No te. — Positio ns of the radio core ( Carilli. Owen fc Harris! 1994 ; Carilli et al 



19971 : ICarillilll995h and the systemic redshifts based on the He II line are adopted 



as the frame of reference in this paper. 



Table 2 

Summary of observations and instrumental setups sorted by position angle (P. A.) 



Object 


Date 


Instrument 


Setup 


Seeing 


Slit 


Kesln. 


A Coverage 


P.A . 


Obs.Time 






(") 


(") 


(A) 


(A) 


(deg) . 


00 


B2 0902+34 


16/01/02 


LRIS 


LS 600/5000 


0.9 


1.5 


4 


5150-7650 


68.7 


2x1800 




07/01/02 


NIRSPEC 


Low Disp. N5 


0.9 


0.76 


14 


14600-17400 


68.5 


6x900 


4C 60.07 


15/01/02 


LRIS 


LS 600/7500 


0.7 


1.5 


4 


5450-7600 


81.0 


2x1800 




25/02/01 


ESI 


Low Disp. 


0.8 


1.0 


13 a 


4000-9600 


135.2 


3.5x1800 




07/01/02 


NIRSPEC 


Low Disp. N6 


0.9 


0.76 


11 


15600-19800 


81.4 


6x900 


4C 41.17 


24/02/01 


LRIS 


LS 600/7500 


0.7 


1.5 


6 


5200-6100 


19.4 


3x1800 




23/02/01 


LRIS 


MOS 300/5000 


0.6 


1.0 


4 


4300-6600 


21.6 


4x1800 




10/12/96 


LRISp 


LS 400/8500 


0.9 


1.0 


8 


5500-9280 


76.5 


28xl200 b 




15/01/02 


LRIS 


MOS 400/8500 


0.7 


1.0 


6 


5300-8000 


81.0 


3x1800 




03/02/97 


LRIS 


LS 600/5000 


0.9 


1.0 


5 


4320-6850 


170.8 


3x1800 




07/01/02 


NIRSPEC 


Low Disp. N6 


0.9 


0.76 


11 


15600-19800 


42.7 


3x900 




07/01/02 


NIRSPEC 


Low Disp. N6 


0.9 


0.76 


14 


15600-19800 


70.2 


5x900 




07/01/02 


NIRSPEC 


Low Disp. N7 


0.9 


0.76 


14 


20300-25000 


42.7 


2x900 




07/01/02 


NIRSPEC 


Low Disp. N7 


0.9 


0.76 


14 


20300-25000 


67.8 


4x900 



a Near the redshifted Lya line at A = 5825 A; the resolution varies roughly linearly from 3 A at 3900 A to 40 A at 11000 A. 
b Spectropolarimetric observations, a detailed analysis of which was presented in lDev et al. I lll997f) . 
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Table 3 

Summary of Emission-Line Measurements 





Object 


Line 


— x — t — 

rest 

(A) 


A obs 
(A) 






-r— 






i-'WHAd 

(A) 


FWHM 

(kms L ) 


(Jontinuum Flux 

,,.-17 

(10 erg s 


Line Plux 

— 1 — 2 \ 
1 cm z ) 


tiw obs 

(A) 


B2 


0902 + 34 (core) 


Lye 


1215.7 


5334.69 ± 0.09 


3 


.3883 


- 


U 


.0001 


21.9 ± 0.2 


1233 ± 13 


18.3 ± 0.5 


7890 ± 95 


431 ± 12 






C IV(doublet) 


1549.0 


6795.00 ± 0.62 


3 


.3807 


± 





.0004 


22.7 ± 1.5 
17.2 ± 0.9 


1003 ± 65 
718 ±"37 


16.3 ± 0.4 


1770 ± 70 


109 ± 5 
127 ± 7 






He If ' 


1640.5 


7199.37 ± 0.37 


3 


.3881, 


± 





0002 


11.8 ± 0.4 


1502 ± 60 


B2 


0902+34 (total) 


Lya 


1215.7 


5336.26 ± 0.23 


3 


.3896 


± 





.0002 


21.4 ± 0.4 


1203 ± 24 


27.4 ± 3.8 


25640 ± 600 


937 ± 133 




C IV(doublet) 
He iY ' 


1549 
1640.5 


6799.58 ± 0.72 


3 


.3897 


± 


o 


.0005 


10.4 ± 1.7 


460 ± 75 


21.6 ± 2.9 


2590 ± 420 


120 ± 25 


4C 


60.07 (core) 


Lya 
C IV 


1215.7 


5826.52 ± 07 


3 


79*>8 


± 


o 


0001 


47 1 ± 2 


2424 ± 12 


< 0.01 


41.1 ± 0.2 


> 5000 






1548.2 


7415.7 ± 0.37 


3 


.7899 


± 





000 1 


81.4 


6x900 












C IV 


1550.8 


7428.1 ± 0.37 


3 


.7899 


± 





000 1 
















Ho II 


1640.5 


7855.7 ± 1.07 


3 


.7887 


± 





.0007 


73.9 ± 2.7 


2820 ± 105 


< 0.01 


9.5 ± 0.3 


> 790 






C III] 


1908.7 


9140.7 ± 1.70 


3 


.7889 


± 





.0007 












4C 


60.07 (total) 




1215.7 


5824.59 ± 0.18 


3 


7913 


± 





.0002 


44.9 ± 0.5 


2310 ± 25 


< 0.01 


56.1 ± 0.7 


> 4000 




c y rv(doubiot) 


1549.0 


7420.75 ± 1.34 


3 


7907 


± 





.0009 


36.5 ±3.3 


1475 ± 130 


< 0.01 


5.0 ± 0.3 


> 600 






He II 


1640.5 


7848.50 ± 4.41 


3 


.7813 


± 





.0027 


84 


3220 ± 4400 


<: o.o2 


4.0 ± 1.6 


5 230 






C III] 


1908.7 


9141.88 ± 1.82 


3 


.7896 


± 





.0010 


59.9 ± 4.7 


1960 ± 155 


^ 0.01 


0.2 ± 0.8 


"> 14 


4C 


60.07 (filament) 


Lya 


1215.7 


5818.07 ± 0.29 


3 


7859 


± 





.0002 


20.0 ± 0.7 


1030 ± 36 


<: o.oi 


6.2 ± 0.2 


> 6600 


4C 


60.07 (nearby galaxy) 


[O III 


3727.0 


6973.30 ± 0.33 





8710 


± 





.0001 


30.0 ± 0.8 


1290 ± 34 


< 0.06 


9.0 ± 0.1 


> 140 


4C 


60.07 (galaxy at filament) 


[O II] 


3727.0 


7046.98 ± 4.25 





.8908 


± 





0011 


51.7 ± ? 


2200 ± 450 


<: o.o2 


3.7 ± 0.2 


!> 180 


4C 


41.17 (core) 


c y rv 


1215.7 


5834.5 ± 0.1 


3 


8002 


± 





.0001 




613/1373 ± 13/45 




146 ± 5.5 








1548.2 


7428.9 ± 0.2 


3 


.7981 


± 





.0001 




541 ± 14 




7.5 ± 0.2 








C IV 


1550.8 


7441.2 ± 0.2 


3 


7981 


± 





.0001 




541 ± 14 




5.7 ± 0.2 








Hell 


1640.5 


7870.7 ± 0.4 


3 


.7979 


± 





0002 




553 ± 28 




5.5 ± 0.3 








C III] 


1908.7 


9152.5 ± 0.6 


3 


7951 


± 





.0003 




511/1120 ± 151/135 




9.1 ± 0.2 








[O II] 


3727.0 


17890.36 ± 0.54 


3 


.8002 


± 





.0001 


42.4 ± 1.1 


711 ± 18 


57.8 ±7.1 










H/3 . 


4861.0 


23351.72 ± 3.27 


3 


.8039 


± 





.0007 


33.0 ± 8.1 


424± 104 


0.19 ± 0.04 


17.2 ± 4.9 


92.7 ± 32.9 



[O III] 4959 23816.29 ± 1.20 3.8026 ± 0.0002 36.1 ± 3.0 454 ± 38 0.18 ± 0.04 59.8 ± 5.2 337.8 ± 81.2 

[O III] 5006.9 24046.74 ± 0.42 3.8027 ± 0.0001 37.9 ± 1.1 472 ± 13 0.18 ± 0.04 222.8 ± 5.4 1260 ± 295 

4C 41.17(total) Kf3 4861 23355.84 ± 8.35 3.8047 ± 0.0017 55 700 ± 270 0.47+0.12 52.5 ± 21.5 112 ± 54 

[O III 4959 23810.62 ± 3.03 3.8015 ± 0.0006 50.8 ± 7.8 640 + 98 0.50 + 0.14 129.4 + 20.5 259 + 85 

[O III] 5006.9 24041.56 ± 1.13 3.8017 ± 0.0002 58.4 ± 3.0 728 ± 37 0.47 + 0.16 484.7+22.6 1031 + 351 



Note. — Where two values are given, line is split. 



